Fiche 201 Niveau : 2 


1 Introduction 


Une des activités principales de l’analyse est l’approximation de fonctions. De ce point de vue le 
théorème le plus célèbre est le théorème d’approximation de Weierstrass : 


Théorème 1.1 L'espace vectoriel des fonctions polynomiales est dense dans l’espace vectoriel C{a, b] 
des fonctions continues sur un intervalle compact |a, b], muni de la topologie de la convergence uni- 
forme sur ce compact. 


Autrement dit, si on considère les fonctions particulières 1,x,x°?,-::,2",--., la fermeture pour la 


topologie de la convergence uniforme sur [a, b] de l’espace vectoriel engendré par ces fonctions est 
l’espace C‘{a, b] tout entier. On dit que la famille (1,x,---,x",---) est totale dans C[a, b]. 

Si on considère maintenant la famille (x,æ,-..,2",...) on peut se poser la question suivante : 
à quelle condition cette famille est-elle totale dans C{a, b] ? Le théorème de Müntz - Szäsz répond 
à cette question. On commencera par établir un théorème de densité pour l’espace L?[a, b] puis on 
établira le résultat pour C{a, b]. 

Il existe diverses démonstrations de ce théorème. Celle que nous proposons ici n’utilise que des outils 
élémentaires (mais elle n’est pas simple pour autant) et suit l’exposé qui est fait sur cette question dans 
le livre de Philip J. Davis intitulé /nterpolation and Approximation (Blaisdell Publishing Compagny 
1963). 


2 Problème 


Partie I - Matrices et déterminants de Gram 


Soit H un espace de Hilbert dont on note (, ) le produit scalaire et || || la norme. Soient x1,%2,--- ,%» 
des éléments de H. On définit la matrice de Gram relative à ce système et son déterminant par : 


(list) (idea) << Cit) 
G(x1,%2,--- nn) — no? És un Fu 
Ci: de) OS 


et 
g(x) = det (G(x, T2, 4) ; 


1) On définit pour 1 < à < n le vecteur : 
nm 
Vi — nn URDEE 
j=1 


1 


on note À la matrice (a;.;);,; et À = (G;;);,; Sa trans-conjuguée. 


Montrer que : 
Gya, V2, ,Yn) = AG(1, T0, + ,Tn) À 
et que : 
g(Y1; Y2, * ** Un) = | det(4)|° g(t1, To, , Ta). 
2) Nous allons montrer quelques propriétés de g(x1,%2,--: ,%à). 
a) Montrer que : 
ne) #0 
et que g(x1,t2,--: ,&n) = 0 si et seulement si le système (11,2, --- ,x,) est lié. 
b) Montrer que : 
g(is a, ++ Ta) € [ml tel? [|sall 
et que l’égalité a lieu si et seulement si le système (21,22, --- ,%,) est orthogonal. 


c) Application à l’inégalité d’'Hadamard. Soit B = (b; ;) ; une matrice à coefficients complexes 
où |b;,;| < M. Montrer que : 
| det(B)| < M'n"/?. 


3) On suppose maintenant que les vecteurs z1,%2,::- ,%, Sont linéairement indépendants. On note 
E = [x:,%2,--- ,x,] le sous-espace de dimension n engendré par ces vecteurs. Si y € H on note 
Ôg(y) la distance de y au sous-espace vectoriel Æ. 

Montrer que : 
(y)? _ g(t1, , "tt; Zn; y) 
g(x1, T2, °°: ) 

et que le vecteur z de E qui réalise ||y — z|| = min,eg |[[y — x|| est donné par : 
En). Arno 265 AR) 0) 
“ (on Di) boy do) Er (aida): (ra, U) 

g PA T2 .... T 
Po à FN San MAD) Mer Mes) 
T] To -.. Tn 0 


Partie II - Le théorème de Müntz - Szâsz 


1) Déterminant de Cauchy. Nous utiliserons dans la suite un déterminant de Cauchy, C’est-à-dire 
un déterminant de la forme : 


1 1 À 

AP a +b2 a1 +bn 

De _ PE a2 se a2 +bn 
mir, poils, las 1 

an+b1 an+b2 an +bn 


où les dénominateurs sont non nuls. Montrer que : 


IL>i(ai — a5)(b — si 


i>j>1l 
D, = J2 


AC + b;) 


2) Soit (pi)i>1 ie suite de nombres réels deux à deux distincts tels que p; > +. Soit q un nombre 
réel tel que q > —> 5. On pose : 


1 

1 : 2 

Ô— min [e-u ax? —...— ax" dr | . 
{ak}g=1ien \Jo 


62 = 1 IT( Pi — q | 
29+125 \m+a+l 


1 


Montrer que : 


3) Soit (x?),e7 une famille infinie de fonctions puissance distinctes où p > —à. Montrer que pour 
que cette famille soit totale dans l’espace L?[0, 1] muni de la norme habituelle 


lie (f 1ropar) 


il faut et 11 suffit que la famille des exposants contienne une suite p; telle que l’une des conditions 
suivantes soit satisfaite : 
: 1 QE 1 
1. lim; p; = —5 et ru (p ni i) = +00. 
2. lim; Pi = p avec —5 < p < +00. 


+00 1 
i=1 pi 


1/2 


3. lim;_,, pi = +00, p; £ Det = +00. 


4) Soit (x?),e7 une famille infinie de fonctions puissance distinctes où p > 0. Montrer que pour 
que cette famille soit totale dans l’espace C[0, 1] muni de la norme uniforme 


FI = ne | f(x) 


x€[0,1 


il suffit que l’un des exposants p soit nul et que la famille des exposants contienne une suite p; telle 
que l’une des conditions suivantes soit satisfaite : 


É +00 
late sine 0er = = +00. 
2. lim;_,, p; = p avec 0 < p < +oo. 


3 Solution 
Partie I - Matrices et déterminants de Gram 


1) Faisons le calcul explicite (et bestial) du coefficient d;; de la matrice À G{x;,%2,-:: ,%,) À. 
Pour cela commençons par calculer le coefficient c; ; du produit À G{x1,%2,-:: ,%n). 


Ge ou Zn Ep) = (5x su) = (Yi, Ti). 
On obtient alors : 
nm 
hs —= > Ga Dr (Us, æ) AO .)- RACE 
S—1 


On obtient bien le coefficient d’indice (4, j) de la matrice G{(y1, yo, : + , Un). 


3 


De cette formule on tire pour le déterminant : 


g(y1, Ye, - HS Un) — det(A)g(x1, x, - Le 4 det(A). 


Mais comme det(À) — det(A) on obtient en définitive : 
g(yi; V2, °°: su) = Idet(A)[° g(æ1; Ta2,°:: RE 


2) Quelques propriétés de g(x1,%2,--: ,æ). 
a) Soit E — [x:1,x2,--- ,x,] le sous-espace vectoriel engendré par les vecteurs 21, %2,--- ,2n. 
Soit FE, un sous espace de dimension n contenant Æ. Considérons une base orthonormée (u1,u2,-:-,u») 
de E,. Alors g(ui,u2,:-: ,u,) = 1. Par ailleurs, tout vecteur x; se décompose sur la base (u;); : 


n 
Ti — ) di,jUj: 
J=1 


Où &;,; = (æ;,u;). Soit À la matrice dont les coefficients sont les a; ;. Cette matrice a pour vecteurs 
colonnes les vecteurs x; (décomposés dans la base (u;);). En conséquence le determinant de À est 
nul si et seulement si les vecteurs x; sont liés. Or : 


gti, T2, °° En) = Idet(A)lg(us, u2,: ++; un) = |det(A)f. 


Donc : 
LL g(t1, &, + Lx) > 0, 
2. g(t1,%2,:-+ , Zn) — 0 si et seulement si les vecteurs x; sont liés. 
b) Soit z1,%2,:-- , x, des vecteurs linéairement indépendants 
(sinon on sait que g(x1,%2,- +: ,%n) — O0). 


Rappelons la méthode d’orthonormalisation de Schmidt : on construit par récurrence : 


= #1 
77 Ten] 
et pour tout j > 2: 
1. E;_; est le sous-espace engendré par x1,%2,--:,%;_ 
(ou ce qui revient au même par y1, Yo, * ** , V1). 


2. z; est le vecteur x; dont on soustrait sa projection orthogonale sur E;-_: : 


25 = Tj — Pr, .(t;), 


puis, 
Uj — . 
= 
11251) 
Ainsi les vecteurs 71, Y2, : :: ,Y, forment un système orthonormé et pour tout 1 < 3j <nona: 
[ei T2,°:: T5] = [un U2,°:: US]. 


De ce fait, la matrice À qui exprime les y; en fonction des x; est une matrice triangulaire inférieure. 
On a donc : 
2 2 , 2 
|det(A)[= la11/1a22l--" lann| 


4 


et donc : 


) : 
T T: .. T EE ——_———_———_—_—_—]— —— — ——_—__ _ —, 
9 1: 2; ) n lai 1/21a22/2 + lan nl? 


Remarquons que le coefficient a; ; est le coefficient de x; dans l’expression de y; en fonction de 
T1,TL2,°°° ,2,. On a donc : 


1 
Qi = © — 
021 
d’où j 
— = ft; — Pr,,(x;)] < [x;l|, 
A5, 


l'égalité ayant lieu, si et seulement si la projection de x; sur E;_; est nulle, c’est-à-dire si et seulement 
si æ; est orthogonal à tous les x; où à < 7. 
On en conclut que : 
2 2 2 
gt, %n) € [a elfe anti 
l'égalité ayant lieu si et seulement si les vecteur x; forment un système orthogonal. 


Remarquons que cette démonstration indique l’interprétation géométrique suivante : le déterminant 
g(t1,to,--: , tn) est le carré de l’hypervolume dans l’espace à n dimensions de l’hyperparallélépi- 
pède construit sur les vecteurs x;. 


c) Soit (ui, U2, - :: , u,) un système orthonormée. On définit les vecteurs 


n 
j=1 


Alors on sait d’après les résultats précédents que : 
OZ gftr, 22, ,2n) = |det(B)F < [il [all [taf 


Mais comme la base (u;); est orthonormée, 


nm 
rl? = D lb, < nM°. 
ÿ=1 


Donc : 
ll real < (RM?) 
et par suite : 
| det(B)| < M'n£. 
3) Soit E = [r1,2%2,--- ,x,]. En vertu de la formule démontrée à partir de la méthode d’orthonor- 


malisation de Schmidt au I-2)b) on peut écrire : 
g(x1, to, on er) — g(t1, T2, "tt; Ln; )Iy _ Pe(y)l. 


Mais ||y — P#(y)||? est le carré de la distance de y au sous espace FE, ce qui prouve la formule : 


(T1, L2,°"" , La, VU 
y? = Lt em U) 
g(ti, To, be) 


Soit z — Px(y) la projection orthogonale de y sur le sous-espace Æ. On a y = z + (y — z) et y — 2 
est orthogonal à Æ. En cnséquence pour tout x; on a 


(y, B;) = (& Ti) 


5 


et ceci caractérise z parmi tous les points de Æ. Si nous écrivons z sous la forme : 


2 = Qi + Q2X2 + **: + Ann, 
alors les coordonnées a, &2, : : - a, Sont la solution du système : 


a1(T1, T1) + Got, To) ++ + An(l1, Ln) = (X1,Y) 


Q(Tn, T1) + Q2(Tn, La) * ** + An(ln, Tn) —= (Ln;Ÿ) 


Les formules de Cramer nous donnent les valeurs des a; sous la forme du quotient de déterminants : 


asér) bite Lee MU) ée tite) 

(to, T1) (T2, To) 7° (T2, y) De (T2, Tn) 
_ lEman) (mag) te (amy) te (em) 
' g(1, T2, --* ,Tn) 


D'un autre côté, si on développe suivant la dernière ligne : 


(T1: Ti) (T1, To) HAE (D) (T1, y) 

Fe. (To, T1) (T2, To) see (Las Ta) (T2, y) 

g(x1, T2,°°° En) ben) A ne a (tn, y) 
Ti To -.. Ts 0 


en tenant compte des valeurs de a; calculées par les formules de Cramer, on trouve bien ÿ 7", a;x:. 
Partie II - Le théorème de Müntz - Szäâsz 


1) Le déterminant de Cauchy. Considérons les a; et les b; comme des variables indépendantes. 
Les dénominateurs sont donc des polynômes homogènes de degré 1, tandis que les numérateurs sont 
des polynômes de degré 0. Le déterminant D), est donc une fraction rationnelle de degré —n. Le 


dénominateur commun est 
nm 


IL («+0 


ij=1 
de degré n?, donc 
Pa(a,--- OA DE bn) 


D; un ; 
IT, (a + b;) 
où P;(@1,-:: , an, 01, ,0,) est un polynôme de degré n°? —n. Si a; = a; (avec à £ j) alors deux co- 
lonnes sont identiques et D,, = 0. Il en est de même si b; = b;. En conséquence P,(a1,--: ,an,b1,--- ,b;) 


est multiple du polynôme : 


n 


LL (@— 0) —b;). 


i>j>1 


Or ce dernier polynôme est exactement de degré n°? — n. En conséquence : 


ESNC — @;)(b: — b;) 
TC + b;) 


6 


Di-—6; 


où c, est une constante. Il reste à montrer que cette constante est 1. Remarquons que : 


1 1 ns 1 
ai+b1 a1+b2 ai +bn 
CE, HET ie 1 
a2+b a2 +b a2 +b. 
an D =. 2 1 2 2; 2 nm 
An An ….. An 
An +b1 An +b2 An +bn 
d’où on tire tout d’abord que : 
1 1 
a1i+b1  a1+b2 a1 +bn 
| +6 +6 + 
lim An D = a2 1 a2 2 a2 n 
An — +00 ; 
1 1 1 
puis : 
1 1 0 
a1+b1 a1 +b2 ai +bn-1 
L 0 
a2+b1 a2+b2 a2+bn-—1 
En: . nt. -@,..D; = ; | 
bn — +00 An — +00 ; 1 
Gn—1+01 an—1+b2 An—1+0n—1 
1 1 1 
Par conséquent : 
, … D Cn 
lim lim il 


bn—+00 an—+00 Dy_] En1 


Comme il est facile de voir que c: = 1 on conclut que pour tout n on a c, = 1. 


2) Les conditions p; > —2 et q > —2 nous permettent de dire que les fonctions x et x? sont 


dans L?[0, 1]. C’est dans cet espace de Hilbert que nous allons nous placer maintenant. Appliquons la 
formule : 


T1,%2, °° ,Ln, 
82 — Op(y)? = g( 1, #2 y) 
g(T1, T2, , Tn) 
Calculons les produits scalaires : 


1 
1 
ne | D 
0 PARA 
En conséquence les deux déterminants g(x1,%2,-:- ,%n, y) et g(x1,t2,--: ,x,) sont des déterminants 


de Cauchy qu’on sait calculer : 


[5-1 (Pi _. p;)* 
TE Pb 


Se EC VE C 
UT M (g+a+1) [lp +p5 +1) [f(g + mm + 1)2° 


g(t1, &, . à — 


et donc : , ; 
6? _ TE (q Fe Di) 
(29 +1)[[521( + pi +1)? 


3) Nous allons étudier ce qu’il se passe dans les trois cas prévus. 


1. lim; p; = —3 et S (pi + À) = +00. Écrivons : 


.+1\2 
[Lg + pi +1)? [l' (1 n piti 


i=1 


Fe : 
Comme p; — —2 on a pour ? suffisamment grand 


et donc le numérateur : 


n 1 
II (: _ Di + 2) 
T 
i=1 g+> 
est borné. D'un autre côté le dénominateur : 


_ pi + à 
nos 


i=1 


diverge vers +00 car la série : 
+00 


pi+s 
D. 


i=1 


est divergente. On conséquence le quotient tend vers 0, ce qui montre que toute fonction x‘ 
(q > 0) est approchée par une suite de points du sous espace engendré par les x?'. Comme les 
7 pour q entier > 0 forment une suite totale dans L?[0, 1], il en est de même de la famille x?:. 


2. lim; Ps = p avec —5 < p < +00. Alors : 


Pi-q _ p—q 


D —© = ———— 
imopitq+tl p+rq+li 


Or pour p > —; et q > Oil est facile de voir que : 


Soit 


Il existe n tel que pour tout à > n on ait: 


Pi — q pP—q ” 
— — — | <eE, 
© DEN Ou 
et donc : 
Pig et 
pitq+i 
Le produit : 


(2a +1) [T(a + pi + 1)? 
converge donc vers Ü ce qui montre que pour tout q > 0 x° peut être approché par une suite 
d’éléments de [x?:]. Comme la famille x? (q entier > 0) est totale dans Z?[0, 1] il en est de 
même de la famille x?:. 


3. lim;_,. p; = +00, p; £ Det pce _ = +00. Dans ce cas on écrit : 


Mut-n%  Mali- sh 
IL(a + pi +1)? TE (i + ei) 


i=1 
Le numérateur reste borné, tandis que le produit du dénominateur diverge vers +co. On conclut 


comme dans les deux cas précédents. 


Réciproquement, si la famille x? est totale on doit montrer qu’on se trouve dans un des trois cas pré- 
cédents. Si la famille des exposants p ne contient aucune sous-suite réalisant l’une des trois conditions 
prévues, c’est que on est dans un des cas suivants : 


1. la famille (p) est une suite p; telle que 
li ee 
Énpere De 
avec 
+00 1 
> ( Ë :) + 
i=1 


2. la famille (p) est une suite p,; telle que 


li = 
RER 
avec 
+00 1 
>» — < +00. 
Pi 


i=1 

3. la famille (p) est une suite p; qui se décompose en deux sous-suites, l’une du premier type, 
l’autre du deuxième type. 

Choisissons un entier q > 0 différent de tous les p;. Les calculs faits dans la partie directe de la 


démonstration montrent que dans aucun de ces cas, x? peut être approché aussi près qu’on veut par 
une combinaison linéaire des xPi. 


4) Soit n > 0 et p; > 0, alors 


1 k Pi— 
<n pe CA UE d 
< Î D « 
On voit par cette dernière inégalité qu’on va faire reposer la bonne approximation dans C[0, 1] des 
éléments de C[0, 1] par des combinaisons linéaires des +,, par la bonne approximation dans L?[0, 1] 
de ces éléments par des combinaisons linéaires des %,,_1. Étudions alors les deux cas indiqués pour 
cette partie : 


— — +00. Avec cette hypothèse les x7-* vérifient la 


condition 3) du résultat de la section précédente. On conclut que la famille x,,,_1 est totale dans 
L?[0, 1]. Conjugué avec l'inégalité précédente ce résultat implique que tout x” où n > 0 peut 
être approché par une combinaison linéaire des x,,,. Si en plus on a un des p qui est nul, alors 
1 est aussi approché (et même atteint), donc en vertu du théorème de Weiersrtrass (les x” où 
n > 0 forment une famille totale dans C0, 1]) les x,, forment aussi une famille totale dans 
C10, 1]. 

. Him; pi; = p avec 0 < p < +0. Dans ce cas si p > 5 on est de nouveau dans un cas 
d'application du résultat de la section précédente. On conclut que la famille x,,_; est totale 
dans Z?[0,1}, puis que les x,, forment une famille totale dans C{0,1]. Si p < + on ne peut 
appliquer directement les résultats de la section précédente car p — 1 < —. On travaille alors 
avec la famille x? où c est une constante > 0 telle que cp; > 5. Cette famille est alors comme 
on vient de voir totale dans C[0, 1]. Comme l’application de [0, 1] dans lui même qui à x fait 
correspondre x° est bijective en faisant le changement de variable x’ — x° on déduit que la 
famille x”: est totale dans C'{0, 1], ce qu’on voulait prouver. 


—1 
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